FUNCTION DecodeEPSG &
!
(epsg) &
!
RESULT (gridmapping)
IMPLICIT NONE
!Arguments with intent in:
INTEGER (KIND = short), INTENT(IN) :: epsg
!Local declarations:
TYPE (CRS) :: gridmapping
REAL (KIND = float) :: centralMeridian
!----------------------end of declarations-------------------------------------
gridmapping % epsg = epsg
SELECT CASE (epsg)
CASE(4326)
CALL SetCRS (GEODETIC, WGS84, gridmapping)
CALL SetGeodeticParameters (gridmapping, prime_meridian = 0.0)
CASE(4230)
CALL SetCRS (GEODETIC, ED50, gridmapping)
CALL SetGeodeticParameters (gridmapping, prime_meridian = 0.0)
CASE(3003) !Gauss Boaga West
CALL SetCRS (TM, ROME40, gridmapping)
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = 9. * degToRad, &
falseE = 1500000., falseN = 0., k = 0.9996)
CASE(3004) !Gauss Boaga East
CALL SetCRS (TM, ROME40, gridmapping)
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = 15. * degToRad, &
falseE = 2520000., falseN = 0., k = 0.9996)
CASE(21781) !CH1903 Swiss topo
CALL SetCRS (SOC, CH1903, gridmapping)
CALL SetSwissParameters &
(gridmapping, latc = 0.819474, lonc = 0.129845, &
azimuth = 1.570796, falseE = 600000., falseN = 200000., k = 1.)
CASE(23028) !ED50 / UTM zone 28N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 28 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23029) !ED50 / UTM zone 29N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 29 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23030) !ED50 / UTM zone 30N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 30 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23031) !ED50 / UTM zone 31N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 31 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23032) !ED50 / UTM zone 32N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 32 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23033) !ED50 / UTM zone 33N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 33 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23034) !ED50 / UTM zone 34N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 34 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23035) !ED50 / UTM zone 35N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 35 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23036) !ED50 / UTM zone 36N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 36 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23037) !ED50 / UTM zone 37N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 37 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(23038) !ED50 / UTM zone 38N
CALL SetCRS (TM, ED50, gridmapping)
centralMeridian = (6 * 38 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32601) !WGS 84 / UTM zone 1N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 1 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32602) !WGS 84 / UTM zone 2N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 2 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32603) !WGS 84 / UTM zone 3N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 3 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32604) !WGS 84 / UTM zone 4N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 4 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32605) !WGS 84 / UTM zone 5N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 5 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32606) !WGS 84 / UTM zone 6N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 6 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32607) !WGS 84 / UTM zone 7N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 7 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32608) !WGS 84 / UTM zone 8N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 8 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32609) !WGS 84 / UTM zone 9N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 9 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32610) !WGS 84 / UTM zone 10N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 10 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32611) !WGS 84 / UTM zone 11N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 11 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32612) !WGS 84 / UTM zone 12N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 12 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32613) !WGS 84 / UTM zone 13N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 13 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32614) !WGS 84 / UTM zone 14N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 14 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32615) !WGS 84 / UTM zone 15N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 15 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32616) !WGS 84 / UTM zone 16N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 16 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32617) !WGS 84 / UTM zone 17N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 17 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32618) !WGS 84 / UTM zone 18N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 18 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32619) !WGS 84 / UTM zone 19N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 19 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32620) !WGS 84 / UTM zone 20N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 20 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32621) !WGS 84 / UTM zone 21N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 21 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32622) !WGS 84 / UTM zone 22N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 22 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32623) !WGS 84 / UTM zone 23N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 23 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32624) !WGS 84 / UTM zone 24N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 24 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32625) !WGS 84 / UTM zone 25N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 25 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32626) !WGS 84 / UTM zone 26N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 26 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32627) !WGS 84 / UTM zone 27N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 27 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32628) !WGS 84 / UTM zone 28N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 28 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32629) !WGS 84 / UTM zone 29N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 29 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32630) !WGS 84 / UTM zone 30N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 30 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32631) !WGS 84 / UTM zone 31N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 31 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32632) !WGS 84 / UTM zone 32N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 32 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32633) !WGS 84 / UTM zone 33N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 33 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32634) !WGS 84 / UTM zone 34N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 34 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32635) !WGS 84 / UTM zone 35N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 35 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32636) !WGS 84 / UTM zone 36N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 36 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32637) !WGS 84 / UTM zone 37N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 37 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32638) !WGS 84 / UTM zone 38N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 38 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32639) !WGS 84 / UTM zone 39N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 39 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32640) !WGS 84 / UTM zone 40N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 40 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32641) !WGS 84 / UTM zone 41N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 41 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32642) !WGS 84 / UTM zone 42N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 42 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32643) !WGS 84 / UTM zone 43N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 43 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32644) !WGS 84 / UTM zone 44N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 44 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32645) !WGS 84 / UTM zone 45N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 45 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32646) !WGS 84 / UTM zone 46N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 46 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32647) !WGS 84 / UTM zone 47N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 47 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32648) !WGS 84 / UTM zone 48N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 48 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32649) !WGS 84 / UTM zone 49N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 49 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32650) !WGS 84 / UTM zone 50N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 50 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32651) !WGS 84 / UTM zone 51N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 51 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32652) !WGS 84 / UTM zone 52N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 52 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32653) !WGS 84 / UTM zone 53N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 53 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32654) !WGS 84 / UTM zone 54N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 54 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32655) !WGS 84 / UTM zone 55N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 55 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32656) !WGS 84 / UTM zone 56N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 56 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32657) !WGS 84 / UTM zone 57N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 57 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32658) !WGS 84 / UTM zone 58N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 58 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32659) !WGS 84 / UTM zone 59N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 59 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32660) !WGS 84 / UTM zone 60N
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 60 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 0., k = 0.9996)
CASE(32701) !WGS 84 / UTM zone 1S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 1 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32702) !WGS 84 / UTM zone 2S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 2 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32703) !WGS 84 / UTM zone 3S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 3 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32704) !WGS 84 / UTM zone 4S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 4 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32705) !WGS 84 / UTM zone 5S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 5 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32706) !WGS 84 / UTM zone 6S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 6 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32707) !WGS 84 / UTM zone 7S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 7 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32708) !WGS 84 / UTM zone 8S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 8 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32709) !WGS 84 / UTM zone 9S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 9 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32710) !WGS 84 / UTM zone 10S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 10 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32711) !WGS 84 / UTM zone 11S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 11 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32712) !WGS 84 / UTM zone 12S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 12 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32713) !WGS 84 / UTM zone 13S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 13 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32714) !WGS 84 / UTM zone 14S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 14 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32715) !WGS 84 / UTM zone 15S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 15 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32716) !WGS 84 / UTM zone 16S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 16 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32717) !WGS 84 / UTM zone 17S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 17 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32718) !WGS 84 / UTM zone 18S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 18 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32719) !WGS 84 / UTM zone 19S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 19 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32720) !WGS 84 / UTM zone 20S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 20 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32721) !WGS 84 / UTM zone 21S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 21 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32722) !WGS 84 / UTM zone 22S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 22 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32723) !WGS 84 / UTM zone 23S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 23 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32724) !WGS 84 / UTM zone 24S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 24 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32725) !WGS 84 / UTM zone 25S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 25 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32726) !WGS 84 / UTM zone 26S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 26 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32727) !WGS 84 / UTM zone 27S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 27 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32728) !WGS 84 / UTM zone 28S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 28 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32729) !WGS 84 / UTM zone 29S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 29 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32730) !WGS 84 / UTM zone 30S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 30 + 177) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32731) !WGS 84 / UTM zone 31S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 31 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32732) !WGS 84 / UTM zone 32S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 32 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32733) !WGS 84 / UTM zone 33S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 33 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32734) !WGS 84 / UTM zone 34S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 34 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32735) !WGS 84 / UTM zone 35S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 35 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32736) !WGS 84 / UTM zone 36S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 36 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32737) !WGS 84 / UTM zone 37S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 37 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32738) !WGS 84 / UTM zone 38S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 38 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32739) !WGS 84 / UTM zone 39S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 39 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32740) !WGS 84 / UTM zone 40S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 40 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32741) !WGS 84 / UTM zone 41S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 41 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32742) !WGS 84 / UTM zone 42S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 42 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32743) !WGS 84 / UTM zone 43S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 43 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32744) !WGS 84 / UTM zone 44S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 44 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32745) !WGS 84 / UTM zone 45S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 45 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32746) !WGS 84 / UTM zone 46S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 46 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32747) !WGS 84 / UTM zone 47S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 47 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32749) !WGS 84 / UTM zone 49S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 49 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32750) !WGS 84 / UTM zone 50S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 50 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32751) !WGS 84 / UTM zone 51S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 51 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32752) !WGS 84 / UTM zone 52S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 52 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32753) !WGS 84 / UTM zone 53S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 53 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32754) !WGS 84 / UTM zone 54S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 54 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32755) !WGS 84 / UTM zone 55S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 55 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32756) !WGS 84 / UTM zone 56S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 56 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32757) !WGS 84 / UTM zone 57S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 57 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32758) !WGS 84 / UTM zone 58S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 58 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32759) !WGS 84 / UTM zone 59S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 59 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
CASE(32760) !WGS 84 / UTM zone 60S
CALL SetCRS (TM, WGS84, gridmapping)
centralMeridian = (6 * 60 - 183) * degToRad
CALL SetTransverseMercatorParameters &
(gridmapping, lat0 = 0., centM = centralMeridian, &
falseE = 500000., falseN = 10000000., k = 0.9996)
END SELECT
RETURN
END FUNCTION DecodeEPSG